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C^ ■ Abstract 



We report sensitive Chandra X-ray non-detections of two unusual, luminous Iron Low-Ionization Broad 
Absorption Line Quasars (FeLoBALs). The observations do detect a non-BAL, wide-binary companion 
quasar to one of the FeLoBAL quasars. We combine X-ray-derived column density lower limits (assuming 
^,^ . solar metallicity) with column densities measured from ultraviolet spectra and CLOUDY photoionization 

\^ ' simulations to explore whether constant-density slabs at broad-line region densities can match the physi- 

os . cal parameters of these two BAL outflows, and find that they cannot. In the "overlapping-trough" object 

(-H ' SDSS J0300-I-0048, we measure the column density of the X-ray absorbing gas to he Nh > L8 x 10^"* cm"^. 

Oj. From the presence of Fe ii UV78 absorption but lack of Fe ii UV195/UV196 absorption, we infer the density 

' ' in that part of the absorbing region to be n^, ^ 10^ cm"^. We do find that a slab of gas at that density might 

^ . be able to explain this object's absorption. In the Fe m-dominant object SDSS J2215-0045, the X-ray ab- 

"^ ' sorbing column density of Nh > 3.4 x 10'"* cm"^ is consistent with the Fe iii-derived Nh > 2 x 10^^ cm"^ 

d . provided the ionization parameter is log L'^ > LO for both the rig = 10" cm"^ and n^ - 10'^ cm"^ scenarios 

considered (such densities are required to produce Fe in absorption without Fe ii absorption). However, the 
p^ . velocity width of the absorption rules out its being concentrated in a single slab at these densities. Instead, 

^ ' this object's spectrum can be explained by a low density, high ionization and high temperature disk wind 

■"nI" , that encounters and ablates higher density, lower ionization Fe iii-emitting clumps. 

■^ ■ Keywords: quasai's: general, absorption fines, individual (SDSS J030000.56-H004828.0, SDSS 

(^ ; J221511.94-004549.9,SDSS J025959.684-004813.6) 

o 

^"T ' 1. Introduction 



Broad Absorption Line (BAL) quasar spectra are not as common as typical quasar spectra but provide 
a unique look at the central regions of Active Galactic Nuclei (AGN). A BAL quasar is characterized by 
absorption troughs fro m gas with blueshifted outflow velocities of typically 10% the speed of light (0.1c) 



d ' (IWeymann et al.lll99 1'). The lower limit on the widt h of a BAL trough is i n part a matter of definition; the 



traditional required minimu m width is 2000 k m s dWevmann et al.lll991l) . but a minimum width of 1000 



km s"' has also been used (iTrump et al.ll2006l) . Such velocity widths are larger than essentially all galactic 
wind outflows, and thus ensure a sample of outflows driven predominantly by AGN; of course, narrower 
AGN-driven outflows can and do exist. 
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BAL quasars themselves are distributed into three subtypes: high-ionization (HiBAL), low-ionization 
(LoBAL), and iron LoBAL (FeLoBAL). HiBAL quasars show absorption from only relatively high-ionization 
species such as C iv, N v and Si iv. LoBAL quasars also show absorption from low-ionization species such 
as Mg II, Al III and Al ii. FeLoBAL quasars are LoBAL quasars with absorption from one or more excited 
states of Fe ii or Fe in. Note that a range of ionization stages is seen even in HiBAL quasars, and that in 
LoBAL and FeLoBAL quasars the ionization simply extends to lower ionization stages JHall et al.ll2002h . 

Approximately 10% of quasars in optically selected samples with spectra covering rest-frame 1400- 
1550 A exhibit a trough or troughs which can be considered broad absorption. The true frac tion is expected to 
be hig her due to selection effects which bias surveys against the detection of BAL quasars. ^ Hewett and Foltz 
(120031) reported a correct ed traditional BAL qu asar fraction of 22 + 4% using their sample of 42 bright 
{Bj < 19) BAL quasars. iReichard et alj (120031) estimated a co rrected tradit i onal BAL quasar fraction of 
15.9 



1 .4% using a sample of 224 BAL quasars with / < 20. iTrump et al.l (12006 ) report an uncorrected 
traditional BAL quasar fraction of 10.4 ± 0.2% using a sample of 1756 BAL quasars with / < 20, or an 
uncorrected fraction of 26.0 ± 0.3% using a less conservative "absorptio n index" criterion t o define a BAL 
quasar (uncertainties from Pois son statist i cs on ly). By comparing the Trump et al.l (120061) catalog to the 
near-infrared 2MASS database. iDai et al.l ( 120081) report a corrected BAL fraction of 43 + 2% for luminous 
quasars (Mk^ < -30.1) using the absorption index criterion. However. iKnigge et alj (120081) find that the 
absorption index BAL quasar criterion includes a large number of false positives. Taking that int o account, 
they find uncorrected and corrected BAL fractions of 13.6% and 17+3%. Scaling the results of Dai et alj 
(|2008) to account for the contamination found by Knigge et al. ( 2008) yields a BAL fraction of 23% , which 
is also the upper limit BAL fraction quoted by Knigge et al. (2008) . Lastly, Ganguly et alj (l2007h find an 
uncorrected BAL fraction of 1 1%, corrected to 23% in Gangulv and BrothertonI (l2008h . In summary, there is 
general agreement that ^23% of quasars exhibit BAL troughs after selection effects are taken into account. 
The fraction of quasars exhibiting BAL troughs could be due to an orientation effect, such that all or 
most quasars have BAL outfl ows covering some of the lines of sight along which they would be seen as 
quasars (e.g.. iTurnshekll 19861) . Alternatively, it could be due to an evo lutionary effect, such that all or most 
quasars have BAL outflows for some fraction of their lifetime (e.g., IS urdei and Hutsemekerslll987b . Of 
course, some combination of the above is also possible (e.g., [Morris 19881). In an y case, these outflows ma y 
have mass loss rates comparable to the accretion rates required to power quasars dSteenbrugge et alj|2005h . 
Thus, understanding BAL quasars and BAL outflows is important for understanding quasars as a whole. 



O bservations show th at BAL quasars are X-ray weak compared with reg ular quasars (e.g. | Green et aL 

19951 [Gallagher et alj|200 6. and references therein). The disk wind model of Murrav et alj (11995! ) explains 
this and other properties of BAL quasars well. In that model, the X-ray weakness is attributed to intrinsic 
absorption from 'hitch-hiking' gas between the BAL outflow and the X-ray emitting region. The nearly 
completely ionized metal atoms in the hitch-hiking gas are kept ionized by capturing electrons and then 
immediately absorbing X-rays from the quasar, just as an H ii region is kept ionized by protons capturing 
electrons and then immediately absorbing photons of energies 13.6 eV or greater Because highly ionized 
metal atoms can only efficiently absorb X-rays, the hitch-hiking gas protects the BAL outflow from overi- 
onization by X-rays while transmitting the ultraviolet (UV) radiation that accelerates the BAL outflow to 
high velocities.' We assume a BAL quasar has an intrinsically normal X-ray spectrum but is shielded by 
some absorbing gas, so that constraints can be placed on the hydrogen column density required to absorb the 
X-ray flux. 

It is also possible that some or all BAL quasars are intrinsically X-ray weak. It has been shown by 



The hitch-hiking gas may also act as an X-ray reflector and a soft X-ray emitter, but we do not consider those effects herein. 



Giustini et alj (120081) that some BAL quasars may be less X-ray absorbed than previously thought, and the 
X-ray weakness is due to a differing spectrum from typical quasars. The authors admit there may be some 
factors that could affect this result. These include the X-ray energy range used for spectral analysis, partial 
covering of the source by the X-ray absorber, or an ionized X-ray absorber. IStreblyanska et al.l (120101) go on 
to show that when an ionized absorber is considered, the column density of intrinsic absorption significantly 
increases. Although the possibility of intrinsic X-ray weakness cannot be fully rejected by current data sets, 
traditional BAL quasars (with BI> 0, where BI is the BALnicity index) appear to be always X-ray absorbed 
dStreblvanska et al.ll2010l) . Both the above studies u se some BAL quasai-s from lTrump et al.l (l2006l) . which 
may not be actual BAL quasars ("Knig ge et al.ll2008h . 

Recent sky surveys have also unc overed a slew of BAL quasars which exhibit unusua l properties as 
compared to typical BAL quasars (e.g.. lHall et al.ll2002tlBrunner et akluOOSHPuc et al.ll2002l) . Investigation 
of these properties is as worthwhile as the study of BAL quasars themselves, as the most unusual objects will 
define the parameter space spanned by BAL outflows in general. If we are to have a complete understanding 
of quasars, our models must explain all normal and unusual quasar behaviour 

In this paper we discuss continued research into two such unusual BAL quasars. SDSS J030000.56-I- 
004828.0 (hereafter SDSS J0300+0048, Fig.[I}i) is an "ov erlapping-trough''BAL quasar at z = 0.89191 with 
nearly complete UV absorption below rest-frame 2800 A (JHaU et al.ll2002l: iHall et aLll2003b . This FeLoBAL 
quasar has a non-BAL binary companion, SDSS J025959.68H-004813.6 (hereafter SDSS J0259+0048), 
which is located 19'.'5 away from SDSS J0300+0048 at a redshift z = 0.894 (Av = 330 + 160km s"'). 

The second object is SDSS J221511.94-004549.9 (hereafter SDSS J2 2 15 -0045 Fig. [II)), a reddened 
FeLoBAL with detached Mgii, Aim, Alii, and Feiii UV34,48 absorption dHall et al.ll2002l) . This quasar is 
atz = 1.4755 as calculated by the associated Mg ii absorption. SDSS J2215-0045 ha s absorption from Fe m 
but there are no absorption feature s from Fe ii, making i t a very rare and unusual find. iHall and HutsemekersI 
(L2003) discuss how the work of .de Kool et al.l (120021) shows that the presence of Fern but lack of Fen 
absorption restricts the outflow to high densities and a narrow range of column densities. 

X-ray data has been obtained using the Chandra X-ray Observatory and analyzed to help understand these 
two unusual BAL quasars. We discuss our data and methods in ^ present the results of our observations 
in ^and JH discus s their implications in ^and summarize our results in ^ We adopt the cosmology of 



Spergeletal.1 (120071) : Ho= 73.2, Qm = 0.259, and Oa = 0.741. 



2. X-ray Data and Analysis 

Observations of our two targets were carried out in Chandra Cycle 4 using the Advanced CCD Imaging 
Spectrometer (ACIS). 

We use the observed-frame energy range 0.5-6.65 keV: 0.5 keV is an obvious lower energy cutoff due 
to the increased background (and higher throughput degradation) at low energies, and 6.65 keV was chosen 
as the upper limit because it has the same effective area as the 0.5 keV lower limit and because the effective 
observing area drops quickly at energies >6.65 keV. The source count rate in the energy range 0.5-6.65 keV 
is expected to be a factor of 100 higher than in the range 6.65-13 keV, while the expected background count 
rate is the same in these energy ranges.^ 

Given the small number of detections of X-ray photons over a finite time period, we are limited by small- 
number statistics rather than by the background level. To determine upper limits at specified confidence 
intervals in the face of small-number statistics, we use the Bayesian method of iKraft et al.i (il991i) rather than 



- See]http://cxc.harvard.edu/proposer/POG/html/index.htnil| 



the frequentist method outlined in lGehrelsl (119861) . The Bayesian method is favored due to its ability to work 



with non-zero background flux when in the zero source count regime. 

2.1. AGN X-ray-UV Luminosity Correlation 

The X-ray and UV luminosities of regular (non-BAL) quasars are observed to have a strong correlation 
across multiple decades of UV luminosity; this offers a useful tool for measuring the X-ray weakness of 
a BAL quasar, which tend to have far weaker X-ray luminosity than a typical quasar. The correlation is 
usually expressed in terms of the logarithm of the ratio between the 2 keV (hkev) and 2500A (/2500) rest- 
frame specific luminosities, denoted as Ogx'- 

aox = 0.372 log(/2kev//2500) (1) 



This correlation has been quantified carefully in a recent study bv lStefFen et al.l (120061) by observing the /; 



1-2 keV 



and ^2500 of 333 AGNs (from multiple surveys). The best fit to the data is: 

iogihkev) = (0.721 + 0.011) log(/25oo) + (4.531 + 0.688) (2) 

Thus, provided a UV luminosity, the exp ected X-ray luminos ity can be calculated. For typical quasars, Oox 
can be written as a function of /2500, as in lSteffen et al.l (120061) : 



OoxiUV) = (-0.137 + 0.008) log(Z250o) + (2.638 ± 0.240) (3) 

To characterize the X-ray weakness we compare the observed ag^ value, calculated using the observed 
^2500 and l2kev in equation ([TJ, to the expected a„x value agxiUV). The difference between the two is denoted 

Affo.v: 

Affo,. = aox - aox(UV) (4) 

By exploiting the strong correlation between Z2500 and l2kev in ordinary quasars, we can quantify the 
X-ray weakness of our targets. 

To determine /250() for our targets requires estimating /25()o for their underlying continua. We used the 
spectra available for each object in SDSS Data Release Six ( Adelman-McCarthvet al. 2008) . We used 
either the smoothed or extrapolated flux at 2500 A rest- frame as /2500; details are given in the individual 
objects' discussions below. We have assumed a value of 10% as the uncertainty on the UV flux for the 
following discussion. The luminosity distance to each quasar was calculated following iPerh (il999i) based on 
the quasar's redshift and our adopted cosmology. 

3. Results: SDSS J0300+0048 and SDSS J0259+0048 

SDSS J0300+0048 and SDSS J0259-i-0048 were observed simultaneously for 6743.5 seconds on Decem- 
ber 19, 2002 (UT). No X-ray photons were detected within a 2'.'5 radius of the position of SDSS J0300+0048. 
SDSS J0259+0048 was detected, with 39 photons within 2'.'5 of its position. The measured background of 
the image is such that we expect only 0.244 + 0.014 observed-frame 0.5-6.65 keV background photons in 
the detection apertures. 



3.1. Results: SDSS J0259+0048 

We first analyze the SDSS J0259+0048 results, as it is a non-BAL quasar and is thus expected to be 
unremarkable. As a test of our methods, we seek to reproduce the observed X-ray photon count of 39, by 
using only the observed UV flux and the relations in ^2.11 

For this quasar at z = 0.894 we find /25oo= (0.55 + 0.05) x 10"^^ erg s"' cm^^ Hz"' and Z2500 = (2.00 + 
0.20) X 10^°ergs-i Hz"'. 

Using equation (|2]i we therefore expect log(/2fev) - -31.2^q7 erg s^'cm"^ Hz"'. Given the predicted 
X-ray flux, we expect 2.8 times the observed count rate of 1.59 in the log, but the difference is < Icr. Thus, 
our method accurately reproduces the observed X-ray photon count of an unremarkable quasar. 

Using equation ([B we calculated an observed Uox - -1.63 + 0.04 for SDSS J0259-I-0048. The expected 
aox is calculated by equation ([3]l. Plugging in the observed /2500 we expect aoxiUV) = -1.51 + 0.34. Thus 
AffoA = -0.12 + 0.34. Therefore, within the uncertainties, the observed and expected values for Cov (and 
X-ray counts) are the same. 

3.2. Results: SDSS J0300+0048 

Using the same method as above, we analyze the observed and predicted X-ray quantities for the 
FeLoBAL quasar SDSS J0300+0048 at z = 0.89194 to determine the column density of the intervening 
matter. 

The "detection " of zero photons in the image results in an upper limit of 3.00 photons at 95% confidence 



dKraft et al.lll991h . equivalent to a count rate of < 4.445 x 10"^ photons s ' in the observed bandpass. This 
corresponds to an observed X-ray flux upper limit of /2fov < 1.793 x 10"-''' erg s"' cm"^ Hz"' at 95% 
confidence. 

To find /25()o we fit a power law in /,, to narrow normalization windows centered at 3060 A and 4735 A 
rest-frame. The resulting flux at 2500 A is f25oo = (3.77 + 0.38) x 10"^^ erg s"' cm"^ Hz"'. 

Equation ^ then predicts the X-ray flux of SDSS J0300H-0048 to be log(f2kev) = -30.58!^:^^ ergs s"' 
cm"^ Hz"', which includes the uncertainty in the power law slope in the UV band. Therefore, we expect to 
see 2.6^Q^ X-ray counts in the log in the image. This is higher than the observed 95% confidence limit of 
0.5 in the log. 

We characterize this X-ray weakness as follows. The expected aox (from equation |3]l is agxiUV) - 
-1.64 + 0.35 and the observed aox (from equation [U is aox ^ -2.35 at 95% confidence. Therefore the 
deviation of the observed aox from the predicted value is Aaox < -0.72, which is a reduction in f2i(ev by 
a factor of 164 from the expected flux. This reduction in X-ray flux is attributed to shielding gas near the 
quasar, which must have a hydrogen column density of Nh > 1 -8 x 10^"* cm"^ at the quasar redshift, assuming 
solar metallicity. That lower limit was calculated assuming neutral gas, but applies to neutral or ionized gas. 
If the gas is partly or fully ionized, it will absorb fewer photons per unit column density, and therefore the 
req uired minimu m column density will be larger than the above lower limit. 

Gibson et al.l (12009:) report aox ^ -2.33 and Aaox ^ -0.70 for this object in their summary of BAL 



quasars from SDSS data release 5. The values are slightly different from values quoted in this paper because 
slightly different formulae were used. Nevertheless, the values are consistent within the uncertainties. 

We can compare our result with the column densities re quired to explain the UV absorption, in particular 
the Ca II column density A^Ca// - (7.13±1.15)xlO''* cm~^ (IHall et al.l2003h . Because its ionization potential 



is less than that of H i, Ca 11 becomes the dominant calcium ion only at large columns behind an H i ionization 



^^http://cxc. harvard.edu/toolkit/pimms.jsp: tlie online version of Portable, Interactive Multi-Mission Simulator, providing count-rate 
estimations and predictions for Chandra. 



front (§5.2 of lHall et alj2003h . iFerland and PerssonI (1 1989k found that for clouds of density 10''^ cm ^ Ca ii 



is dominant only at a column densities Nh > 6.3 x 10^"* cm"^, more than a factor of ten higher than the 
column density of the H i ionization front. 

We have used the photoionization simulator CLOUDY^ to investigate the absorption in various ions for a 
representati ve Broad Line Regio n (BLR) (density of rie = 10^' cm"^, log L'^ = -1.5, where U is the ionization 
parameter) JBaldwin et alj2003b . Note that our CLOUDY simulations were run with solar metallicity and no 



dust obscuration. Even at Nh = 10^^ cm~^, Ca ii is not the dominant calcium ion. However, the ionic fraction 
of Ca II immediately behind the H I ionization front is ~\%forng - 10" cm"^, significantly higher than the 
~0.3% seen for«p = 10^'^ cm"^. Given that the abundance of Can is -5.64 in the log, each A^;/ = 10^^ cm"^ 
behind the H i ionization front yields 10(^2-2-5.64) y^orth of Ca ii: Ncaii - 2.29 x 10'"* cm"^. Ca i has a similar 
behavior to Ca ii though its ionic abundance immediately behind the H i ionization front is only 0.0001 %. 

Thus, to match the Can column observed in SDSS J0300+0048 with standard BLR parameters requires 
only Nh = 3.11 x 10^^ cm"^ behind the hydrogen ionization front. This agrees with the constraint on the 
total Nh behind the front from the observed upper limit on the Ca i column density: Nh < 7 x lO^"* cm"^. 

However, the UV-derived Nh - 3.11 x 10^^ cm"^ is much less than the X-ray-derived Nh > 1.8 x 
lO'^'* cm"^. This can be explained by assuming a higher ionization parameter at the cloud face. Using 
CLOUDY, we investigated the behavior of Ca ii with higher values for the ionization parameter; namely 
log U - 0.5,1.0,1.5, and 2.0. We also varied the density: separate runs for a density of tie - 10" cm "^ and 
n,, = 10'^ cm"^ were carried out for each of the above values of log U. The results are shown in Figure|2] 

The shaded regions are defined by the UV-derived Ca ii column density (horizontal region) and the X-ray 
derived lower limit on the hydrogen column density (vertical regions). The lower limit on the hydrogen 
column is shown at Icr (lightest grey), 2cr, and 3cr (darkest grey). The curves represent the different log U 
of interest for this object. For a specific log U to match the observations, its curve must fall in the allowed 
region for both the column densities plotted. In Figure |2l the allowed region for each log U is represented 
by a cyan highlight. Thus, this object requires log U > 0.0 for both n,, = 10" cm"^ and n^ - 10'^ cm"^ 
(Figures |2^ and|2]3, respectively). This lower limit is consistent with the constraints posed by Cai (plotted 
in Figure O, and also the constraints from absorption in Mg i, Mg ii, and Fe ii. 

However, as we shall see in § 15.11 observations of Fen* rule out the log {/ > 0, «g - 10" cm"^ scenario; 
gas with such parameters may form part of the absorber in this object, but it cannot form the entire absorber. 

4. Results: SDSS J2215-0045 

SDSS J2215-0045 was observed for 6615.9 seconds on June 21, 2003 (UT). The measured background 
of the image is such that we expect only 0. 102 + 0.004 observed-frame 0.5-6.65 keV background photons in 
the detection apertures. 

To find /2500 we fit a power law in fy to narrow normalization windows centered at 2668 A and 3035 A 
rest-frame. The result is /250() = (7.26 ± 0.73) x 10"^^ erg s"' cm"^ Hz"'. The true /2500 could be somewhat 
higher since we do not account for possible Mg 11 absorption in the shorter- wa velength normalization region , 
nor do we account for the inferred continuum reddening of £'(B-y)^0.06 cHall and HutsemekersI 120031) . 



Accounting for either of those effects would increase /2500, which would increase the predicted /2fey, which 
would in turn increase the absorption needed to bring the predicted f2kev below the upper hmit we measure. 
Thus, this value of /2500 is conservative with respect to the amount of X-ray absorption we infer 



http://www.ferland.org/cloudy/ 



No photons are detected in the observed 0.5-6.65 keV bandpass, so we have a conservative upper Hmit 
of 3.00 photons at 95% confidence.^ This limit corresponds to a count rate of < 4.544 x 10""* photons s"' . 
The observed X-ray flux limit is therefore f2kev< 1 .776 x 10"^^ erg s"' cm"^ Hz"' . From the observed X-ray 
flux upper limit we calculate the observed aox upper limit, using equation ([U, and find Oox < -2.45 at 95% 
confidence. 

The expected X-ray flux is log(/2fev) = -30.53;^^^^ erg s"' cm"^ Hz"', which includes the uncertainty 
in the power law slope in the UV band. Given that flux, 2.70^^ g^ counts in the log should be observed over 
6615.9s, which is higher than the observed 95% confidence limit of 0.5 in the log. We also expect ao^tiUV) - 
-1.74 + 0.35. The deviation of the observed q-qi from the predicted value is Aq-qi < -0.71. Compared to 
SDSS J0300+0048, we probe higher-energy photons in SDSS J2215-0045 due to its higher redshift, and so 
to produce the same observed Aaox < -0.71 requires a higher absorbing column: Nh > 3.4 x 10^"* cm"^ at 
the quasar redshift, assu ming solar metallicity. Again, this is a hard lower limit assuming a neutral absorber 

Gibson et al.l (120091) also evaluated values for this object. T hey find a„r < -2.1 \ and Ao-ov < -0.96. 



Again, due to slightly different formulae, the values derived in iGibson et al.l (I2009h are slightly different 
than those calculated here. In this case, the Auox value is < 1.5cr different. 

The presence of doubly ioni zed iron (Fe iii) absorp tion without any significant singly ionized iron (Fe ii) 



absorption is a rare occurrence. Ide Kool et al.l ( 120021) showed that this iron absorption behaviour could be 



produced with a density of n^. = 10 cm"^, an ionization parameter log U - -1.5, and a very narrow 
range of total hydrogen column density around log Nh ~ 22 A. We investigated this result with CLOUDY 
simulations using the same parameters as for SDSS J0300-I-0048. The results for Fe in are located in Figure 
|4]and the results for Fe ii in Figure|5] 

In Figure m the shaded region is shaped by the lower limits of both the X-ray derived hydrogen column 
(vertical lower limits) and the UV-derived Fern column density (horizontal lower limits; details given in 
§ I5.2| |. The curves must fall inside this region to correctly describe the observations. The curves must also 
satisfy the upper limit on the Fen column density (also discussed in § 15. 2b shown in Figure |5] The shaded 
region, in this case, is defined by an upper limit Fe ii column. 

The cyan segments in Figure |4]represent the interval over which each curves match the observations by 
falling within the shaded regions of both Figures |4] and |5] We can place a lower limit of log U > 1 .0 on the 
ionization parameter for both the rie - 10" cm"^ and n^ - 10'^ cm"^ cases (Figures|5t and|5j5, respectively). 
This lower limit is consistent with constraints imposed by Mg n, Al ii, and Al m. 

The physical basis behind these constraints is that Fe m absorption without significant Fe n absorption is 
seen when the absorber has just barely enough column density to form a hydrogen ionization front, because 
Fe m is only abundant just before such a front and Fe ii is abundant just after one. Furthermore, to create an 
Fe m column matching the observations in the narrow regio n just before the front requi res a high density, 
because such high densities force Fe iv to recombine to Fe iii (IBautista and Pradhanlll998l) . 



5. Discussion 

We have shown that both SDSS J0300+0048 and SDSS J2215-0045 are undetected in snapshot X-ray 
observations. Despite the short exposure times of the X-ray observations, non-detections of these optically 
bright quasars require very large absorbing columns and small levels of unabsorbed X-ray scattering (< 1 .3% 
at 95% confidence).^ Note that our column density lower limits assume neutral absorption. As shown by 



'One photon with energy 12.7 keVis detected within 2'.'5 of the position of SDSS J221 5-0045, but photons with such high energies 
are much more hkely to be background photons, as discussed in 



Streblvanska et alj (120101) . when an ionized absorber model is used, BAL quasar X-ray absorbing columns 



can be 1-2 orders of magnitude larger than in the neutral case. Thus, the column densities we quote here are 
hard lower limits, and the true column densities are almost certainly larger. We now discuss the implications 
of these results in more detail. 

5.1. SDSS J0300+0048 

In SDSS J0300+0048, at outflow velocities of 2000 < v < 4000 km s ' we see Mg ii, Mg i, Ca ii and Fe ii 
absorption. At 4000 < v < 10850 kms ', we see Mgn, Fen and excited-state Feii (Feii*) absorption from 
gas which must be located closer to the quasar than the Ca ii-absorbing region. (If it were located farther 
fro m the quasar than the Ca ii-absorbing region, it woul d be shielded by that gas and would itself show Ca ii.) 

Hall et alJ (l2003h and lHaU and Hutsemekerl (|2004 suggested tiiat the SDSS J0300-H0048 outflow could 



be produced by gas in an accretion disk wind if we are looking across the wind rather than down the wind. 
In this model, high-ionization gas is launched from closer to the quasar than low-ionization gas is, with Ca ii 
and Mg i being launched from outside a hydrogen ionization front. ^ 

We can constrain the density in the Feii*-absorbing region by studying the critical densities of excited 
levels from which we either see or do not see Fe ii* absorption. We see absorption from the Fe ii UV78 
multiplet, whose lower level has a critical density of ^ 10^^ cm""*. We do not see absorption from Fen 
UV195/UV196 multiplets, with relevant critical densities of ^ 10^^ cm"^. Thus, we can constrain the 
density in the Fe n*-absorbing region to be ^ jq^"'^*'-'' cm"-'. 

We now examine the ionization behaviour of a slab of gas at this density, to see if a single-slab model can 
explain both the X-ray and UV absorption. In Figure|6]we present two CLOUDY runs for rig = 10^ cm"^, 
with log f/ = 1.0 and logU = -0.5. For log f/ - 1.0, the observed Fen column (Figure |6^) can be 
reached with reasonable hydrogen column densities, but the Can column (Figure |6]3) can only be achieved 
with Nh ~ 10^^ cm"^. When the hydrogen column approaches ~ 10^^ cm"^ between us and the quasar 
UV emission region, the electron scattering optical depth rises above t^s ~ 2, increasing the already high 
intrinsic optical and UV luminosity of this object. (Of course, the X-ray absorption may cover only the X-ray 
continuum-emitting region, and not the larger UV continuum-emitting region, but at the moment we are only 
considering whether a single uniform absorber can explain the data.) For log U = -0.5, the observed Can 
column (FigurelSJ)) is reached at ~ 3 x 10^'* cm"^, at which point the predicted Fe n column is well above 
the observationally inferred lower limit and the temperature within a factor of two of the observationally 
inferred upper limit. 

Thus, in light of the X-ray data, we find that we can nearly match the properties of both the X-ray and UV 
absorbers in SDSS J0300+0048 with one slab of gas of constant density tig - 10^ cm"^ with log U = -0.5 
at its ionized face, implying an absorber at a distance of ~60 pc from the black hole with a thickness of 
~ 3 X 10"* cm or ~1 pc 

If the absorbing gas in SDSS J0300+0048 can be approximated by gas with a uniform density of n^ = 
10^ cm"^, it is difficult to reconcile the resulting ~60 pc distance of the gas from the black hole with the 



^These X-ray scattering constraints are broadly consistent with optical polarization measurements for both objects (Hall, Smith, et 
al., unpubUshed): SDSS J0300+0048 has R-band P = 1.58 ± 0.10% and SDSS J2215-0045 has a white-ligh t P = .40 ± 0.06%. 

^The radial dependence of the density in the outflow at launch is unclear As mentioned in Hall et all J2003h . in region (a) of a 
IShakura and Sunvae\l (1973) disk the density increases with radius. A disk wind launched from that region would have a higher density 
in the Can region than in the Fe ii* region. A low temperature (T < 1 100 K) would be required in the Ca ii region to avoid exciting Fe ii* 
there. On the other hand, the low-ionization absorption in this object covers the Mg ii (and Fe ii) broad emission line region, located 
at radii where the Shakura and Sunyaev ( 1973) disk model predicts a radially decreasing density. If the disk wind is launched at these 
larger radii and its density is decreasing with distance, the density in the Ca ii-absorbing region should be the lowest in the outflow. A 
density of n^ < 10^ cm"^ would be required to avoid populating the lower levels of Fe ii* transitions which are not observed. 



scenario to a disk wind discussed in Hall et al 



original disk wind scenario of iHall et alj (12003 ). Such a distance is more consistent with the alternative 



(120031) : namely, that the outflow sweeps up gas at large 



distances from its origin and decelerates in the process. In this swept-up-gas scenario, the Fe ii*-absorbing 
gas would be gas which has been compressed to n^ > lO-' cm"-' and accelerated to 4000 < v < 10850 km s"' 
by a high-ionization, low-density wind (seen at velocities up to at least v = 10850 kms"' in this object). 
The Ca ii-absorbing gas would be gas which has been overtaken by the wind more recently and which has 
therefore been compressed and accelerated to a lesser degree than the Fen*-absorbing gas. However, it 
remains to be seen (through detailed outflow modeling beyond the scope of this investigation) whether gas 
can be accelerated to 2000 < v < 4000 kms"' without being com pressed to a density «g > lO-' cm""* 
or heated to a temperature T > 1100 K. Meanwhile, as suggested bv lHall et al.l (12003). this swept-up-gas 



scenario can be tested by looking for a long-term increase of the Ca ii BAL outflow velocity. (No increase 
was seen over a rest-frame time span of up to 205 days in multiple SDSS spectra plus our VLT spectrum of 
this object.) 

5.2. SDSS J22 15 -0045 

The conservative upper limit of l ogA^peu < 14.7 (Figure |5ll f or this object was estimated using the 
plausible unabsorbed continuum from lHall and Hutsemekersll2003l (hereafter HH03), shown in Figure lb. 
The details of the estimation of the lower limits for the Fe m column density shown in Figured are given in 
the Appendix. 

We have found in ^ that if a single absorber is responsible for both the observed UV and X-ray ab- 
sorption in SDSS J2215-0045, it must have Nh > 3.4 x lO^"* cm~^ logn^ > 9.5 cm^^ and log U >\, with 
full coverage of the X-ray source and at least 50% coverage of the UV continuum source. We now consider 
whether an absorption system with those properties is in fact plausible. 

If we take Nh = 3 x 10^^ cm"^ and logn^ = 9.5 cm"-', then the absorber is only lO''' cm thick. A larger 
Nh would make for a thicker absorber but would also increase the electron scattering optical depth above 
Tes - 2, incre asing the already hig h intrinsic luminosity of this quasar. The X-ray continuu m region is ~10'^ 



cm in radius (^ Chartas et al.ll2009h . while the UV continuum region is ~10'^ cm in radius (Kochaneket alj 
|2007). To cover both regions, the absorber would have to be at least 20 times larger in the transverse 
direction than it is in the line-of-sight direction. Moreover, the line-of-sight velocity width of the absorber 
is 12,000 km s"\ If that is a turbulent Av, the absorber would double its thickness and halve its density on 
a timescale equal to the crossing time of ~10 days, but observations of SDSS J2215-0045 over rest-frame 
timescales much longer than that reveal no changes in its absorption troughs (HH03). If that Av is a coherent 
velocity spread along our line of sight, the acceleration required to produce it in 10'^ cm is ~100 m s"^. For 
comparison, in the disk wind model of Murray et al. (1995) the radiative acceleration is <100 cm s"^. 

It is hard to see how to bring the above extreme parameters into the realm of the plausible. For example, 
a supersolar Fe abundance would reduce the Nh required to match the observed A^Peiii, and Lya pumping of 
Fe m UV34 (Johansson et al. 2000) would reduce the amount of Fe needed to explain a given A^Peiii- Those 
effects could eliminate the Tes problem, but would shrink the line-of-sight width of the absorber and make 
the dynamical problems worse. 

Thus, a single, uniform absorber with the properties required to explain the UV and X-ray absorption 
in this system is extremely unlikely to exist. Instead, it is likely that the X-ray absorber is Compton-thick 
(> 3 X 10^"^ cm"^) but compact (radius ~ 10'^ cm), so that at most a small part of the UV continuum region 
is obscured by it. The UV absorber must be larger (radius ~ 10'^ cm) and have a column density just under 
that required to form a hydrogen ionization front (Nh % U x \Q-^ cm"^) so that it contains very httle Fen. 
We can constrain -1 < log U < I for the UV absorber The upper limit comes from setting the electron 
scattering optical depth to the UV continuum source to be at most unity. The lower limit comes via our 



measured log A^peiii > 16.62 and logA^peii < 14.7; given those values. Table 2 of Ide Kool et al.l (120021) and 
our Figures |5]and|4]show that log Ug > 9.5 cm"^ and log U > -\^ 

We can reconcile the high density required for the Fe in column with the large velocity and velocity 
spread of Fe iii if we assume the Fe in absorption arises in dense clum ps embedd e d in a lower density, higher 



ionization wind seen in Civ and Siiv, along the lines suggested by IVoit et al.l (119931) . The clump widths 
along the line of sight would be li < NHlne, or d < ip'^ 5±i o (.p^^ T^i& max imum size is larger than the 



typical size of 10'^ cm predicted for putative BLR clouds dKorista et al.lll997b . To cover half the r ~ 10'* 
cm UV source, a large number of such clumps (~ 10^^^ if the clumps are roughly spherical) are required. 

In this object, Fe in absorption is seen at 6000 < Vpeiu < 18000 km s ' along our line of sight, whereas 
Civ and Siiv absorption is seen at < v < 25000 + 3000 kms ' (Figure 1 of HH03). Given that the 
troughs in this object have higher outflow velocities and velocity widths than the average BAL trough, it is 
reasonable to assume that much of the acceleration of the gas occurs along our line of sight; in other words, 
the outflow streamlines are largely parallel to our line of sight. If we imagine an accelerating, relatively low 
density wind (seen in C iv and Si iv) colliding with dense Fe iii clumps, the wind will form a shock around 
the clumps and can in principle accelerate, ablate, compress and bypass the clumps. The acceleration may 
help explain the Fe in trough velocity, along with radiative and possibly magnetohydrodynamic acceleration; 
the ablation helps explains the Fe in trough width; the compression helps explain the inferred high densities; 
and the bypassing ex plains the C iv an d Si iv trough extending to both smaller and larger outflow velocities 
than the Fe in trough jVoit et al.iri993b . 



In summary, our picture of this quasar is that the X-ray-emitting region is completely covered by an X- 
ray absorber and that the more distant UV-emitting region is nearly or completely covered by a UV-absorbing 
wind seen in C iv and Si iv absorption. This wind is accelerated (at least in part along our line of sight) from 
km s"' to 25000 + 3000 km s"' . At some distance or range of distances from the quasar, the wind includes 
dense, Fe in-emitting clumps (with overall covering factor ~50% and low volume filling factor) which have 
been accelerated by the wind to velocities of 6000 < Vpein < 18000 kms"'. A paraboloidal shock surrounds 
each clump, where the wind gas encounters the slower-moving clump gas. The shock-heated gas is no longer 
visible in C iv or Si iv, but between a clump and its shock lies a region where gas ablated from the clump 
will absorb in those ions, at velocities v ^ Vpein- The wind downstream from the clumps — consisting of 
shocked gas and, if the volume filling factor of the clumps is small, unshocked gas as well — is not seen to 
recombine to observable ionization stages when shadowed by a clump, so its recombination timescale must 
be longer than the time it spends in a shadow. 

In this picture, the Fe in-emitting clumps are being ablated and so the Fe in trough depth will decrease 
with time unless the wind encounters new clumps. Such clumps should first appear at low velocity; from their 
absence, we conclude that the wind is not currently encountering new clumps. Spectroscopic monitoring of 
the quasar to look for weakening Fe in troughs (as well as any changes in Fe in trough velocity or velocity 
width) would be worthwhile. 

Lastly, we note the possibility that this quasar may have developed from a typical FeLoBAL. In such 
an object. Fen is seen at low velocities and Civ at both high and low velocities, consistent with a high- 
ionization wind which has encountered clouds optically thick in the Lyman limit and has only begun to 
accelerate them. As time progresses, the clouds will be accelerated and ablated. Eventually, they may be 
optically thin in the Lyman Umit (and thus seen in Fe in and not Fe n) and will be found at higher velocities. 



'^We assume that the UV absorber is fully shielded by a dust-free X-ray absorber. In that case, the X-ray absorber will absorb a 
relatively small fraction of hydrogen-ionizing photons, and the structure of the absorber's Stromgren sphere will be much less aifected 
than the structure of its partially ionized zone. We can therefore use our constant-density CLOUDY simulations to infer the conditions 
of the Fe iii clumps, even though the higher-ionization regions in the outflow must have a lower density. 
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6. Conclusion 

SDSS J0300+0048 and SDSS J2215-0045 have some of the highest known column densities in their 
BAL outflows. 

SDSS J0300+0048 specificafly has a very high Mgii absorbing column, which creates the overlapping 
absorption below rest-frame 2800A. Exploiting the correlation between X-ray and UV luminosities, we find 
h.aox < -0.72, indicating strong X-ray absorption. The X-ray absorbing gas must have a column density of 
Nh > 1.8 X 10^"* cm"^ (assuming solar metallicity). Absorption from FeiiUV78 but not FeiiUV195/UV 196 
in the UV spectra indicates a density in this absorbing region of n^ ^ 10'' cm"^. We have provided two 
CLOUDY runs at this density with ionization parameters log f/ = 1.0 and log U - -0.5. Most properties 
of the absorber in SDSS J0300-I-0048 can be matched with a slab of gas of constant density tie - 10^ cm"-* 
with log U - -0.5 at its ionized face (implying a distance of ~60 pc from the black hole) and a thickness of 
3 ~ 10i^cm(~lpc). 

For SDSS J22 15-0045, we showed that a single absorbing slab of density log «(. > 9.5 cm"-' and log U > 
1 that fully covers the X-ray source and covers at least 50% of the UV continuum could explain the total 
hydrogen column and Fe in-Fe iii column densities. These parameters are shown to have little physical 
plausibility, as they require an abnormally high acceleration to produce the observed line-of-sight velocity 
spread seen in the UV absorption troughs (the velocity spread would otherwise disperse the slab and reduce 
its density on extremely short timescales). One possible scenario to explain the absorption would be dense 
Fe Ill-emitting clumps inside a lower density UV-absorbing wind (seen in C iv and Si iv absorption). The 
UV wind is accelerated from km s"' to 25000+3000 km s"', and impacts the the denser Feiii clumps. 
The clumps are accelerated to the observed velocity dispersion while being ablated by the UV wind. Such a 
model could be tested by looking for evaporation of the dense Fe iii absorbing clumps manifesting itself as 
time-varying absorption features in the UV spectra. Explaining the rarity of Fe m-dominant BAL quasars is 
easier: only a small paramet er space for strong Fe m absorption with little Fe ii absorption present is found 



by both lde Kool et all ( l2002h and our own CLOUDY runs. 
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sity of Wyoming Department of Physics and Astronomy and the Aspen Center for Physics, where parts of 
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NSERC; PBH and MSB acknowledge financial support from Chandra (proposal 04700849). This research 
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cm"' A" . a. The overlapping-ti'ough object SDSS J0300+0048. Below 2800 A, Fe ii and other troughs overlap to absorb nearly all the 
flux from the quasar, b. The Fe m-dominant BAL qua sar SDSS J2215-0045 (solid) and a partial estimate of its intrinsic, unabsorbed 
spectrum (dashed), from lHall and Hutsemekera J2003h . The UV34 and UV48 multiplets of Feni are very strong in absorption, while 
only upper limits can be put on absorption from any Fe n multiplet. 
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Figure 2: The behaviour of the Ca ii column density for varying ionization parameter and density. Shaded regions represent constraints 
posed by column densities of Can and hydrogen in SDSS J0300+0048. The vertical lines separating the dark, light, and lightest 
shades of grey are the lower hmits for the X-ray calculated hydrogen column density for the 90%, 95% and 99% confidence intervals, 
respectively. The hght grey horizontal bar represents Ca ii column density and its uncertainty. The log of the ionization parameter is 
labeled beside its corresponding curve. The cyan(black in hardcopy) highlights represent the regions for each ionization parameter 



where they are in the allowed parameter space, (a) Graph for iig 
figure is available in the online version] 
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Figure 3: These curves correspond directly to tlie curves seen in Figure |2] as a clieck of how consistent the lower limit on the log of 
the ionization parameter is with other ions. A hard upper limit (calculated via UV spectroscopy) is set for the Cai column density in 
SDSS J0300+0048. The same hydrogen column density is used for the vertically shaded regions. Again, the curves are labeled with 
their corresponding ionization parameter, (a) Graph for n^ = lO' ' cm"', (b) Graph for rig = lO'^ cm"'. 
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Figure 4: Curves represent the behaviour of the Fe iii column density for varying ionization parameter and density. The lower limit on 
the hydrogen column in SDSS J2215-0045 is shown at Icr (lightest grey) and 2cr (darker grey). An uncorrected (lightest grey) and 
corrected (darker grey) lower limit on the Fern column is presented (see j|5.2t . The log of the ionization parameter is labeled beside its 
corresponding curve. The cyan (black in hardcopy) highlights denote the regions for which the curves satisfy both the lower limit Fe ni 
column density and the lower limit hydrogen column density as well as the upper limit Fen column density (see Figure|5). Red(light 
grey in hai'dcopy) highlights only satisfy the constraints posed by the Fein and Fen column densities, disregarding our lower limit on 
the hydrogen column, (a) Graph for n^ = 10" cm"'', (b) Graph for n^ = 10'^ cm^^^. [A colour version of this figure is available in the 
online version] 
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Figure 5 : The relationship of Fe m and Fe n absorbing gas in this object is vital. Here we plot the same curves for each ionization 
parameter as seen in Figure |4] but in comparison to the upper limit on Fen absorbing column (see j) l5.2t . The lower limit on the 
hydrogen column density is also the same as seen in Figure|4] These plots show that the ionization parameters supported by the lower 
limit on the Fern column are also supported by the upper limit on the Fen column density, (a) Graph for «(. = lO" cm"'', (b) Graph 
for«, = lO'^cm"'. 
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Figure 6: Two CLOUDY runs for log U = -0.5, 1.0 and rig = 10 cm"^ . The curves representing each ionization parameter are labeled. 
The shaded areas represent the allowed column densities of (a) Fe ii and (b) Ca ii in SDSS J0300+0048 based on UV observed column 
density limits and the lower limit, X-ray derived hydrogen column. 
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Appendix A. Feiii Column Density Lower Limits for SDSS J2215-0045 

The lower limit for the Fe in column density in SDSS J2215-0045 was calculated (after normahzing by the unabsorbed continuum 
shown in Figure lb) using absorption from Fe iii ions in the lower levels of the UV34 and UV48 triplets (3.73 and 5.08 eV above ground, 
respectively). The individual components of each tiiplet are severely blended. Therefore, we treat each triplet's trough as coming from 
a single hne with a Af value equal to the average of the components of the triplet, and then divide the resulting column density by 
three. We assume that the covering factor of the line is unity, so that the optical depth Tv at each velocity v is given by r„ = - log /„, 
where /^ is the absorbed continuum at v normalized by the assumed unabsorbed continuum at v. The resulting r,. is always optically 
thin (Tv < 0.6). We find a raw lower limit of log WfcIii > 15.18 (lightest grey region in Figure|4). Even this is the larges t Fein column 
observed in any outflow to date and is nearly ten times the Np^m observed in FBQS 0840+3633 bv lde Kool et alj j2002l) . 

We can set an improved log A'pein lower limit by attempting to account for Fe in ions in other levels. Assuming unsaturated troughs 
and that the excited Fe in level populations follow a Boltzmann distribution, the relative number densities in the two observed states 
yield Tj, = 34, 250 K. Accounting for Fe in ions in all levels with excitation energy < 5.08 eV at that temperature yields a factor of 28 
correction from observed to total column in Fe in. Thus the lower limit we adopt is log Wpclll > 16.62 (darker grey region in Figure|4). 

It is important to note that the true Fe in colum n density could be substantially higher than any of the above estimates if the troughs 
are saturated, which is quite likely jHall and Hutse mekers 2003). For example, M. Bautista (personal communication) finds a maximum 
possible ratio of A'uv48/A'uv34 = 0.158 at «j > lO'^ cm"^^, whereas we observe a value of 0.447 assuming unsaturated troughs. If we 
conservatively assume that only the UV48 trough is unsaturated and use the maximum possible value of A'uv48/A'fcI1i = 0.00155, again 
found at n, > 10 cm" (M. Bautista, personal communication), we obtain logA'pciii > 17.49. Note that a density n^ = 10 cm"^ 
would imply logA'pcin = 18.5 and a density n,, = 10*-'' cm"' would imply logA'pciil = 19.5. 
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